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Abstract — Solar air heater performance can be enhanced by 
adding roughness to the inner periphery. The present study is 
based on the various shape parameters for broken arc shape 
roughness element of heat transfer and friction factor 
characteristics of the rectangular duct roughened with repeated 
square cross-section broken ribs with a gap. Duct has width to 
height ratio (W/H) of 5, relative roughness pitch (P/e) of 8, 
relative roughness height (e/D) of 0.045 and the arc angle (a) of 
60°. In the range of 0.5-1 gap width (g/e) and 0.16-0.66 gap 
position (d/w) are varied. Under the similar flow condition heat 
transfer and friction factor of smooth duct is compared with 
roughness duct. The duct has Reynolds number (Re) range of 
3000-22300. The maximum enhancement in Nusselt number and 
friction factor is obtained by 6.74 and 6.73 times. 


Index Terms — Heat Transfer, Friction factor, Broken arc, 
Solar energy, Reynolds number, Nusselt number. 

I. Introduction 

The convective heat transfer coefficient between air and the 
absorber plate is low. Convective heat transfer is low due to 
presence of sub-layer and the thermal performance of solar air 
heater is poor. Which can be broken by artificial roughness on 
heat transfer surface [1]. Conversion of solar energy into 
thermal energy is the easiest way for heating application with 
the use of solar air heater. By artificial roughness friction 
losses lead to excess power requirement for the fluid to flow 
through the duct. The turbulence must be created only in a 
region very close to the heat transfer surface to break the 
viscous sub-layer for augmenting the heat transfer. The effect 
of rib shape angle of attack pitch to rib height ratio (P/e) on 
friction factor and heat transfer characteristics for rectangular 
shape with two side roughness wall.. 

Lots of study on heat transfer has been performed on 
several years [3-5]. An idea about some roughness in heat 
transfer can be lead to turbulence flow like situation in the 
region and issue of low heat transfer coefficient can be 
addressed reasonable. Solar radiation incident is absorbed by 
solar collector and it converted into heat for heating the fluid 
(air or water). Their low efficiency is negative effect which is 
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reasonable due to low convective heat transfer coefficient 
between air flow and absorber plate by which increased 
absorber plate temperature. Yadav et al.[7] studied the 
influence of Reynolds number on the surface having array 
externally shape on heat transfer coefficient distribution. The 
enhancement in Nusselt number and friction factor is 
absorbed to be 2.5 times of that of smooth surface. Sethi et 
al.[8] has carry one experimental investigation on dimple 
shape roughness used in solar air heater. The relation 
developed between Nusselt number and friction factor Kumar 
et al.[9] carried experimental research of friction and heat 
transfer in the flowing air in rectangular duct having v-shape 
ribs with gap roughness on the absorber plate. By 
experimentally Nusselt number and friction factor increased 
by maximum 6.74 and 6.73 times. Here correlation also 
developed between Nusselt number and friction factor. Karwa 
and Chitoshiya[10] has experimental carry study of 
performance of thermo hydraulic of solar air heater with 60° 
v-down discrete rib roughness of the absorber plate on the 
side of air flow along with smooth duct solar air heater. The 
enhancement of Nusselt number was maximum2.3 times and 
characteristics of friction were 2.83 times compared with 
smooth duct. Hans et al[12] carried out the effect of multiple 
v-ribs over the absorber plate and generate the relation 
between Nusselt number and friction factor. 

The investigation carried out the heat transfer rate and 
friction factor is purpose to heating of air in rectangular duct 
having multiple broken arc shape as roughness geometry. By 
experimental the relation developed between Nusselt number 
and friction factor. The investigation is given in figure 1. By 
patching aluminium wires in multiple arc shape fashion the 
roughness geometry has been created under the absorber 
plate. Bywhich heat transfer characteristics of rough duct 
have been investigated 

The investigations carried out the heat transfer rate 
and friction factor is purpose to heating of air in rectangular 
duct having multiple broken arc shaped as a roughness 
geometry. By experimental the relation developed between 
Nusselt number and friction factor. The recently investigation 
is given in figure 1. By patching aluminium wires in multiple 
arc-shape fashion the roughness geometry has been created 
under the absorber plate. By which heat transfer 
characteristics of rough duct have been investigated, 25 
roughened plates were fabricated and tested. The obtained 
values are expressed in terms of roughness pitch (p/e), relative 
roughness height (e/D), arc angle (a) and roughness width 
(W/w). All the obtained values are shown in table 1. 
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Nomenclature 

Ap surface area of absorber plate (m 2 ) 

Ao area of orifice meter (m 2 ) 

Cp specific heat of air (J/Kg-K) 

Cd coefficient of discharge 

D hydraulic diameter (m) 

e rib height (m) 

b/e relative broken width 

e/D relative roughness height 

f friction factor for roughened duct 
fs friction factor for smooth duct 

h height of duct 

K thermal conductivity (W/m-k) 

L length of test section in duct (m) 

Nu Nusselt number 

Nus Nusselt number for smooth duct 

APo pressure drop across the orifice plate (pa) 

APD pressure drop across the test section (Pa) 

P pitch (m) 

p/e relative roughness pitch 

Pr Prandtl number 

Qu useful heat gain 

Re Reynolds number 

To air outlet temperature (K) 

Ti air inlet temperature (K) 

Tpm average temperature of absorber plate (K) 

Tfm average temperature of air (K) 

V mean flow velocity in duct (m/s) 

W width of duct (m) 

W/w relative roughness width 
W/H aspect ratio of duct 
Greek symbols: 
a arc angle (°) 
p density (kg/m 3 ) 

P ratio of orifice diameter to pipe diameter 


II. EXPERIMENTAL DETAILS 

2.1. Materials and Procedures 

A rectangular channel have forced convection flow is used in 
this type of experimental set up. The entrance, test and exit 
section are 800 mm, 1000 mm and 500 mm respectively and 
total length of rectangular channel is 2300 mm[l-2]. Wooden 
rectangular duct, electric heater, GI pipe, control valves, 
blower, orifice plate, U-tube manometer, micro -manometer, 
variable transformer, ammeter, voltmeter, thermo couples and 
milli-voltmeter are the components of this experimental setup 
shown in fig.2. The dimensions of inner section are 2300 mm 
x 330 mm x 30 mm. Having cross-section of 1000 mm x 330 
mm an electric heater is provided. Heating wire on a 5 mm 
thick asbestos sheet is obtained by arrangement in series and 
parallel loops. For uniformly relation between absorber plate 
and electric heater thickness of 1 mm mica sheet electric 
heater wire is provided. Variation of heat flux in between 0 to 
1000 W/m 2 is used for variable transformer. A glass wood is 
kept inside the 6 mm thick wooden panel for decrease top loss 
from heater as insulating materials, 12 thermo couple for plate 
temperature measurement and 8 thermo couple for inside 
temperature measurement of duct. For the measurement of the 
temperature total 20 thermo couples are provided. U-tube 
manometer is connected to orifice meter for measurement of 


mass flow rate of air in the rectangular duct. A valve is used to 
control the flow, is called control valves. A digital voltmeter 
is provided for output measurement of the thermo couples. 
Micro manometer is also used to measure pressure drop in test 
section and its least count is 0.01 mm of water. Micro 
manometer, U-tube manometer, milling- voltmeter, voltmeter 
and ammeter are used for finding the function the function is 
in proper way or not and check to all joint, there is no any 
leakage in the system. After verified, blower switched on. By 
control valve of Reynolds number, mass flow rate is arranged. 
Reading was taken in quasi steady process. Quasi- steady 
means slow changes in system that means temperature will be 
steady for 10-15 minute. In order to Reynolds number 
different values of mass flow rate of air was obtained. The 
following parameters are fond by experiment: pressure drop 
across orifice plate, pressure drop across test section, inlet air 
temperature, outlet air temperature and the temperature of the 
plate. 
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Fig.l. roughened multiple broken absorber plate. 

The equation are used for finding the values of heat transfer 
coefficient 4 h’, useful heat gain 4 Q U , Nusselt number 4 N U ’, 
Reynolds number 4 Re’, friction factor 4 f and thermo 
hydraulic performance parameter 4 rf: 

Mass flow rate m has been found by pressure drop across the 
orifice plate 


( 1 ) 


A/r ^ a l a P jPo 

“ %l 1 - 

where coefficient of discharge C d is 0.605 by calibrate. 

And AP 0 = 9.81p m Ah o sin0. Where 0 = 90° for this 

( 2 ) 

Hence, 

AP 0 = 9.81p m Ah 0 (3) 

Friction factor is calculated by 
f = 2p wg Ah d /4pLV 2 

f=2AP D D/4pLV 2 (4) 

Where AP D = p w gAh d , Ah d is head loss across the test 

section length of lm. 

And 
D = ■ 


2{W+ff\ 

Nusselt number (Nu) is calculated by 
Nu = hD/k (5) 

H = Q u /A p (T pm -T fm ) (6) 

where, 4 Q U ’ is rate of heat given by air, 4 A p ’ is the area of the 
absorber plate, 4 T fm ’ is average value of the air temperature 
and 4 T m ’ is absorber plate temperature. 


Q u =mC p (T 0 -T i ) (7) 

From above all these calculation, air properties corresponding 
to bulk mean air temperature is used. By Kline and 
McClintock[14] method value will be calculated. By the 
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analysis of error estimation through various instruments, the 
values of non dimensional number are calculated: Reynolds 
number ±2.20%, Nusselt number ±4.90%, friction factor 
±3.82% and Stanton number is ±3.67%. 



I 



Section at X - X 


Fig2. Schematic diagram of experimental set up. 
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3. Exit Section 
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transformer 

6. Orifice meter 

7. U-tube manometer 


EH-Electric heater H-Duct 
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8. Control valve 
9. Blower 
lO.Electric motor 
11. Selector switch 
12. Variable 

13. Milli-voltmeter. 


2.2. Validity Test 

Compare the value of Nusselt number and friction factor by 
experimental data for smooth duct to the obtained value from 
Dittus-Boelter[15] and modified Blasius equation[16]. 
Nusselt number for smooth rectangular duct is obtained by 
Dittus-Bolter equation as 
Nu s =0.023flp^Pr IJ ' 4 

(8) Friction factor for smooth duct is obtained by modified 
Blasius equation 

f s = 0.085/Ts _<tIS 

(9) 

Comparison of the experimental and estimated 
value for Nusselt number as a function of Reynolds number as 
shown in fig. 3 (a) and friction factor is function of the 
Reynolds number in fig. 3(b). The average variation in Nusselt 
number from experimental value is obtained from equation(8) 


and average variation in friction factor from experimental 
value is obtained from(9). The accuracy of data is obtained by 
experimental set up of these two sets of values. 



Fig.3. (a) Experimental and estimated values of Nusselt 
number of smooth duct. 



Fig.3.(b)Experimental and estimated values of friction factor 
of smooth duct. 


III. Result and discussions 

The heat transfer and friction factor of rectangular duct 
having a heated wall and roughness with multiple broken arc 
shape rib is found by experiments for various roughness 
working parameters. Effect of various parameters on Nusselt 
number and friction factor comes in this section. 

3.1. Reynolds number 

Figure. 3 (a) and 3(b) shows the effect of Reynolds number 
(Re) and friction factor (f) on the Nusselt number (Nu). For 
multiple broken arc shape roughness geometry have pitch 
roughness (P/e) = 8 and other parameters are relative 
roughness height (e/D) = 0.045, arc angle (a) = 60° and 
relative broken width (W/w) = 1 kept constant. These result 
also compared with smooth duct result. Figure 3(a) shows the 
value of Nusselt number is increases in all Reynolds number. 
Due to this experimental value with the increase in Reynolds 
number, turbulence increases which leads to increase the heat 
transfer whereas friction factor decrease with increase in 
Reynolds number, which is shown in figure 3(b). Suppression 
of laminar sub layer for fully developed flow in the duct, 
pressure drop along the duct increases which leads to high 
pumping power. 
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The comparison experimental values and 
predicted values of Nusselt number and friction factor is 
shown in fig.4. This graph shows the values for smooth plate 
in a single graph with lines. 



Fig.4. Comparison of experimental and predicted values of 
Nusselt number and friction factor for smooth plate. 


Fig5. Shows the value of Nusselt number as a function of 
broken width (g/e) for different relative broken position (b/w) 
for 60° inclined ribs duct roughness at selected Reynolds 
number. By observing at any relative broken width the 
Nusselt number is the height for relative broken position of 
0.25 for all Reynolds number. 


P/e=8,g/e=0.5,a=60°,e/D=0.045 
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Fig-5. Effect of relative broken width and relative broken 
position on Nusselt number at fixes Reynolds number. 


In fig(6) and fig(7) shows the Nusselt number and 
friction factor with respect to relative roughness height(e/D) 
for different range Re from 2200 to 22300. The other 
parameters are W/w = 5,a=60°,p/e=8. From the conclude that 
the value of Nu and f increases with increase in e/D for all the 
values of Reynolds number (Re) due to increase in e/D. The 
fact is relative roughness height value increases, roughness 
geometry produces more into the flow which causes more 
turbulence result is increase in Nusselt number (Nu) and 
friction factor(f). The diagram also shows that the rate of 
increase of Nusselt number is lower than that of friction factor 
(f). The fact is higher value of e/D the reattachment free shear 
might not occur and the rate of heat transfer (Q) is not 
proportional to frictional factor. 
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Fig6. Effect of Nusselt number and e/D at various Reynolds 
number. 



Fig.7. Effect of friction factor and e/D at various Reynolds 
number. 


compared the Nusselt number of with or without broken 
multiple arc of the inclined rib with that of the smooth duct. 
The value of Nusselt number ratio(Nu/Nu s ) for different value 
of relative broken position shown in fig(8) and relative broken 
width is shown in fig(9). The ratio of Nusselt number 
increases with increase in relative broken width up to 1.0 
beyond which it decreases with increase in relative broken 
width. Nusselt number ratio is high for relative broken width 
1.0 and low for 2.0, because increase the relative broken 
width beyond 1.0, the flow velocity through the broken will 
reduces, which not be strong to accelerate the flow through 
the broken and then the heat transfer due to this flow not be 
increase continuously in compared with that of the continuous 
ribs. When reduces the broken width lower than 1.0 (i.e b/e = 
0.5). It may leave very small space for flow of the fluid 
through it, which is low turbulence and hence reduce the 
enhancement of heat transfer. 


p/e=8,a=60° 
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Fig.8. Effect of relative gap position on Nusselt number ratio. 
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Fig.9. Effect of relative gap position on Nusselt number. 


The effect of relative broken position (b/w) on the 
Nusselt number ratio for a fixed relative broken width (b/e) of 
1.0 is shown in fig. For any Reynolds number it can be 
obtained the Nusselt ratio is higher for a broken continuous 
rib as compared with without broken rib and the Nusslet 
number ratio increase with increase in relative broken 
position from 0.25 to l.o, attains a maximum at broken 
position of 1.0 and after this it decrease with increase in 
relative broken position. The Nusselt number ratio lies in 
between 1.7-2.58, under similar condition. The variation of 
Nusselt number with relative gap position it presented in fig. 
By which brings out the effect of gap position. Phenomena of 
flow is in rectangular duct a continuous inclined rib gives rise 
to secondary flow along with rib length, which gives the 
working fluid to travel from leading edge to trailing edge of 
the rib. Flow along the ribs is heated continuously and the 
boundary layer grows thicker. The flow turns downwards 
from the side wall and it completes the recirculation loop. 
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Fig.10. Effect of friction factor and W/w at Various Reynolds 
number. 

Fig(10) shows the relative roughness width(W/w) on 
friction factor(f). The friction factor have been plotted is 
function of W/w for various Reynolds number (Re). The other 
parameters are e/D=0.045, a=60°and p/e=8. It has been 
observed that friction factor (f) increases with increase in 
relative roughness width (W/w) and its maximum value of 
W/w value of 7. The angling of wire help in formation of 
secondary flow which promotes turbulence mixing and then 
Nu increases. A single arc along width of duct which having 
W/w value of 1 and considerable improve in heat transfer of 


duct as a compared with smooth one. By increase the number 
of arc would increases the number of secondary flow leads to 
increase in heat transfer and friction factor. The value of W/w 
beyond 5 may result in separation of flow from ribs surface 
and generate secondary layer which reduces heat transfer. The 
friction factor increasing on basis of formation of vortices due 
to separation of flow. 

Fig(ll). Further investigated mechanism of 
enhancement of heat transfer as a result of creating gap, 
ingestion also had been carried out on transverse rib 
roughness surface with or without a gap of ribs. For the study 
take two relative gap position (g/e)=0.25 and 0.5 with relative 
gap width 1.0, have compared in heat transfer for the case of 
continuous rib. Fig(10) shows the variety of Nu with Re for 
transverse ribs with or without gap. In this fig seen that there 
is no significant change in Nu for continuous rib with a gap. It 
shows that the creation of gap in transverse ribs does not 
result significant increase in heat transfer in case of transverse 
ribs. Increase in heat transfer in case of inclined ribs can be 
fully attributed to present of secondary flow. However Falu et 
al[ ] reported that a 90° discrete rib arranged in staggered 
manner on the rough plane with a gap between the ribs 
element by which heat transfer is 10-15% compared to 90° 
continuous rib arrangement. The reason is enhancement of 
heat transfer due to the separation of fluid through the end of 
the discrete rib. 
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Fig.ll. Effect of Nusselt number at Various Reynolds 
number. 


3.2. Friction factor: 

Axial flow profile is distributed by the secondary 
flow measurable effect by which increase the friction 
coefficient in non-circular duct [1]. The lower value of 
friction factor has been repeated for discrete ribs compared to 
the continuous ribs [18], due to abbreviate secondary flow 
cells. Effect of the relative broken width on the friction factor 
of rough ducts with Reynolds number at broken position of 
0.25 is represent in fig [12] in the term of friction factor ratio 
(f/f s ). By which cleared that the friction factor ratio increases 
with increase in Reynolds number due to increase in the 
turbulence flow, while this ratio increase with increase in 
relative broken width upto 1.0 and decrease with further 
increase in the relative broken width up to 2.0. For continuous 
inclined ribs lies close d to friction factor ratio for broken 
width 2.0, due to very weak flow through this large gap (b/e = 
2). The maximum value of the Nusselt number was obtained 
at broken width of 1.0 with a relative broken position 0.25. 
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Fig.12. Effect of broken width on friction factor ratio with 
Reynolds number. 


Correlation for Nusselt number and Friction factor: 

By experimental data it has been observed that Nusselt 
number and friction factor is function of system and other 
parameters are Reynolds number (Re), relative roughness 
pitch (p/e), relative roughness width(W/w) and relative 
roughness height(e/D). 

Correlation for Nusselt number 

The Nusselt number and Reynolds number have power 
relationship is shown in fig. From analysis it is found that 
average slope of all power lines is 1.31. This power line 
equation represented as 

Nu = A 0 ( Re) 1 * 31 

( 10 ) 

The coefficient of A 0 be the function of other influencing 
parameters. Now parameter e/D taking into consideration and 

Nil 

the value of 1 31 =A 0 with respect to e/D graph will be 
plotted. 


( 11 ) 

For second order curative is obtained 

Mm 

Log ( f ^i = logC 0 +C 1 (log(p/e))+C 2 (log(p/e)) 2 

This equation rearrange in 

^,^ 0,3 = Co ( $*** [exp (-0.7 1 (ln( J}} 2 )] 

( 12 ) 



Fig.14. Plot of Nu/tf^ 1 31 with relative roughness height (e/D). 
Where, C 0 is the function of <J>/10. Here the normalizing angle 
of 10° represents the point where maximum in Nusselt 
number. 


Nil 


f _ R<l i.i! i ^2.36 [exp J 


= C 0 


The equation can be rearrange by 

Nil 

( *■ 3 1 ( *■* 3 2 - “ [«*P (— *71. (in (*))*)] 



Fig.13. Plot of Nusselt number (Nu) vs Reynolds 
number(Re)for range of relative roughness pitch. 


Do( A)-"“[exp(-l.S0(ln(i)) i ‘)] 

(13) 

The value of coefficient is 

A 0 =7.46 x 10 -4 , B 0 =3.9 x 10 -3 , C 0 =1.9 x 10" 4 , D 0 =1.89 

x 10 ” 4 

These result in following relationship for the Nusselt number 
are 

Nu = 1.89 x 

10 - 4 (7!s 0 L31 y D - 43 f^r E6 [expC-O.Tian)^} 2 )] (p - M1 ® [ex 
pt-lAOtin^)) 2 }]. (14) 

Fig. 10. Plot of predicted value vs. Experimental value of 
Nusselt number. 
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Fig.15. Plot of predicted value vs. Experimental value of 
Nusselt number. 


Correlation for Friction factor 

Similarly, developed a correlation for the friction factor. 
An=0.15, Bn=4.11, Cn=12.57 

And the final correlation for friction factor can be written as 
f = 12.44 (/Iff} -0 ' 1 " 7 ( -}“ 6 ( £)-«* ( -f A7 

(15) 



Fig. 16. Plot of friction factor(f) vs Reynolds number(Re). 



(e/D) 

Fig.17. Plot of f/Rs ” D 1£ vs relative roughness height (e/D). 



(p/e) 

Fig.18. Plot of t7 Re~ DM (e/Dy sg vs relative roughness 
pitch(p/e). 



( 0 / 10 ) 

Fig.19. Plot of f/Re "°' 19 (ff vs ($/10). 



Fig.20. Plot of predicted value vs. Experimental value of 
friction factor. 
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IV. CONCLUSIONS 

By experimental observation on 60° inclined rectangular duct 
rib roughened with and without broken can conclude that: 

1 . The increase in Nusselt number and friction factor is 

in the range of 1.47-2.57 times and 2.25-2.9 times of 
the smooth duct, respectively for the range of 
Reynolds number 3,000 to 22,300. 

2. The maximum enhancement of Nusselt number (Nu) 

takes place at Reynolds number (Re) of 22,300, 
relative roughness width at 5, relative roughness 
height at0.045, relative roughness pitch 8. 

3. The maximum friction factor obtained at Reynolds 

number of 22300, relative roughness width of 7, 
relative roughness height of 0.045 and relative 
roughness pitch of 8. 

4. Correlation for Nusselt number and friction factor 

have been function of rib spacing, rib height, rib 
angle and Reynolds number. From these correlation 
found that predicted values within the error limit of 
+17% and +12%. 
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